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Abstract 

Recently the Hamilton- Jacobi formulation for first order constrained systems has 
been developed. In such formalism the equations of motion are written as total 
differential equations in many variables. We generalize the Hamilton- Jacobi for- 
mulation for singular systems with second order Lagrangians and apply this new 
formulation to Podolsky electrodynamics, comparing with the results obtained 
through Dirac's method. 

PACS 04.20.Fy- Canonical formalism, Lagrangians and variational principles. 
PACS ll.10.Ef - Lagrangian and Hamiltonian approach. 
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1 Introduction 

Systems with higher order Lagrangians have been studied with increasing interest 
because they appear in many relevant physical problems. As examples we have 
the consistent regularization of ultraviolet divergences in gauge-invariant super- 
symmetric theories Q or effective Lagrangians in gauge theories . Besides this, 
the fact that gauge theories have singular Lagrangians is in itself a motivation to 
the study of the formalism for second order singular Lagrangians. 

The Lagrangian formulation for constrained systems can be found in references 
[Q] and PJ while the Hamiltonian formulation of singular systems is usually made 
through a formalism developed by Dirac |5|, || |7j. In this formalism the constraints 
caused by the Hessian matrix singularity are added to the canonical Hamiltonian 
and then the consistency conditions are worked out, being possible to eliminate 
some degrees of freedom of the system. Dirac also showed that the gauge freedom 
is caused by the presence of first class constraints. 

The study of new formalisms for singular systems may provide new tools to in- 
vestigate these systems. In classical dynamics, different formalisms (Lagrangian, 
Hamiltonian, Hamilton- Jacobi) provide different approaches to the problems, 
each formalism having advantages and disadvantages in the study of some fea- 
tures of the systems and being equivalent among themselves. In the same way, 
different formalisms provide different views of the features of singular systems, 
which justify the interest in their study. 

Here we generalize the Hamilton- Jacobi formalism that was recently developed 
|@, Q to include singular second order Lagrangians. We start in Sect. 2 with an 
overview of the features of singular systems with second order Lagrangians and of 



the Dirac's Hamiltonian formalism for them. In Sect. 3 we develop the Hamilton- 
Jacobi formulation for a general second order system and apply this formalism to 
the case of a singular second order system in Sect. 4. An example is solved using 
both Dirac's and Hamilton- Jacobi formalism in Sect. 5, while Sect. 6 is devoted 
to the conclusions. 

2 Singular systems with second order Lagran- 
gians 

The treatment for theories with higher order derivatives has been first developed 



by Ostrogradski [H| and allows to write the Euler-Lagrange equations, introduce 
conjugated momenta and develop a Hamiltonian formalism for such systems. 
Here we center our attention on the case of a system described by a Lagrangian 
containing time derivatives of the coordinates up to second order (a second order 
Lagrangian) . 

In such a theory the configuration space is described by the n generalized coor- 
dinates qi and its first and second derivatives with respect to the time parameter 
t (we consider a discrete system to simplify the calculations; the generalization 
to continuous systems is straightforward and will be done in Sect. 5). 

The Euler-Lagrange equations, which are obtained from the action integral 

(1) S = f L(q,q,q,t)dt 



using the Hamilton's principle, are given by: 

U dq~ ~dt [dJJ + [djj ~ 

We construct the phase space by introducing the generalized momenta 

dL d ( dL\ 



(3) Pl 



dki dt\dqi 



(4) 7Tv = — 

(conjugated respectively to qi and 9j) and writing the accelerations 9j as functions 
of the coordinates q, velocities q and of the momenta p and tc (q\= fa (g^, q\,pi, 7Tj)). 
The phase space will then be spanned by the canonical variables (qi,Pi) , 7Tjj 
where 9i=<7i. 

By introducing the canonical Hamiltonian defined as 



(5) # c = Pi Qi +Ki q, 



9i=fi L \ q i=fi 



we can write the equations of motion of any function g of the canonical variables 
as: 

(6) 9={g,H c } 

But this procedure is only possible if the determinant of the Hessian matrix 

8 2 L 



(7) H { 



({, i) q, 



does not vanish, otherwise it will not be possible to express all the accelerations 
q\ as functions of the canonical variables and there will be relations as 

(8) $ a (quPi;^,^) = 0; a = l,...,m < 2(n - 1) 



connecting the momenta variables. As a consequence we will not be able to treat 
the canonical variables as an independent set and we have to retreat to formalisms 
specially developed to deal with the dependence among the canonical variables, 
i.e. a formalism for constrained systems. 

The usual treatment of singular systems was developed by Dirac || §, 0] 
who, in order to deal with this problem, introduced a generalized Hamiltonian 
formalism (details can be found in references || and ]TT|]). Such formalism 
can also be applied to the case of second order Lagrangians as can be seen in 



references [j Lj and [1^]. Dirac's formalism consist in considering constraints given 



by equation (|j) as weaA; equations, called primary constraints, and represented 
as: 

(9) $ a (qi,Pi]9i,TT t ) « 

By weak equations we mean those that can't be used until all Poisson brackets 
have been calculated. 

We may add to the canonical Hamiltonian He any linear combination of 
the primary constraints and define a new Hamiltonian, called total Hamiltonian, 
given by 

(10) H T = H c + u a <5> a , 

where arbitrary coefficients. Physically He and H? are equivalent and 

we cannot distinguish between them. The equation of motion for any function 
/ Qi, Pi j 71 i) is given in terms of H T as: 

(11) /- {/, H T } = {/, H c } + u a {/, $ Q } 

The constraints will produce consistency conditions because they must be 
valid at any time and consequently their time derivative must be weakly zero. 
The consistency conditions are given by: 

(12) {$^, H T } = H c } + u a $ Q } w 0; a, (3 = 1, m 

These conditions may be either identically satisfied (when we use the primary 
constraints), determine some of the arbitrary coefficients u, or generate new con- 
straints that will be called secondary constraints. The constraints that have 
null Poisson brackets with all other constraints are called first class constraints 
otherwise they are called second class ones. This classification is completely in- 
dependent of the division in primary and secondary constraints. The extended 



Hamiltonian is defined as 

(13) H E = H C + Vx*x 

were the ty\ include all first class constraints. V\ are arbitrary coefficients and 
we use ([13]) instead of the total Hamiltonian fllQ|) in the equations of motion. 



3 Hamilton-Jacobi formalism for second order 
Lagrangians 

Recently a new formalism for singular first order systems was developed by Giiler 
[§, § who obtained a set of Hamilton-Jacobi partial differential equations for 
such systems using Caratheodory's equivalent Lagrangians method and wrote 
the equations of motion as total differential equations. 

In this section we will use Caratheodory's method to develop the Hamilton- 
Jacobi formalism to a general second order Lagrangian. This formalism can be 
applied to any second order Lagrangian and is not limited to singular ones. The 
singular case will be considered in the next section. 

Caratheodory's equivalent Lagrangians method to second order Lagrangians 
says that, given a Lagrangian L(ft, ?j, t), we can obtain a completely equivalent 
one by: 

/ ■ \ dS(qi,qi,t) 

(14) L — L (ft, q h ft, t) K — >- 

These Lagrangians are equivalent because the action integral given by them 
have simultaneous extremes. So we can choose the function S(qi,Qi,t) in such a 
way that U becomes an extreme and then we reduce the variational problem of 
finding extreme for the Lagrangian L to a problem of differential calculus. To do 
this we must find a set of functions <^(ft, Qi,t) , (ft, t) and S^ft, Qi, t) such that 

(15) L'(ft,A,^,*) = 



and for all neighborhood of k%— A (ft, t) and q\— y?i(ft, ft, t): 

(16) L'(ft,ft,ft,*) >0 

With these conditions satisfied the Lagrangian L' will have a minimum in 
q\— Pi (ft,*) and ft= ipi fft,ft,t) and consequently the action integral will have 
a minimum. So, the solutions of the differential equations will correspond to 
extremes of the action integral. 

From the definition of V we have: 

dS(q h q h t) dS(q i ,q i ,t) dqi dS (ft,ft,t) d ft 



(17) L' = L ft, ft, ft,* 



dt 



dqi 



dt 



dq l 



dt 



;is) 



Using condition ( |i5|) we have: 

dS (ft,ft,i) 



L (ft, ft, ft 



dt 

dSfquQut) . dS(q l ,q l ,t) .. ' 
Qi Qi 



dqi 



dqi 



(19) 



dS_ 
~dt 



<ii=<Pi 



, dS(qi,qi,t) . dS(q h qi,t 
L(q h qi,q h t) ^ '- qi ^— '- ft 



dqi 



dqi 



Since ft= and ft= tfi are minimum points of V we must have 



(20) 

(21) 

or 
(22) 



du 



d ft 



<ii=<Pi 



dL d fdS s 



d q f d qi \ dt 



'dL 




.3 ft 


dQi\ 



dS 



0. 



dqi 



dL 



dft 



ii=Pi 



(23) 



Analogously we must have 
dL' 



dq l 







dL d (dS s 



dq { dqi\dt 



Qi=(3i 



0, 
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(24) 

(25) 

or 
(26) 



dL d dS dS 



dt d ki dqi 
dL 



dS_ 



d 2 S . 

: ?i 

d dS 



d 2 S 



dS_ 



1i=j3i 



dtdq, 
dL d dS 



dqtd <lj 
0, 



<ii=<Pi 



9i=A 



dq^ dtd <li 

From these results, using the definitions for the conjugated momenta given 
by equations @ and (|) and writing Qi—Qi, we have from equation (|19"D that, to 
obtain an extreme of the action, we must get a function S(qi,Qi,t) such that: 



(27) 

(28) 
(29) 



dS_ 
~dt 



H 

Pi 



Pi Qi +VTj q { -L 

_ OS 
dq^ 1 dq { 



These are the fundamental equations of the equivalent Lagrangian method, 
equation ( |2?| ) being called the Hamilton- Jacobi partial differential equation, or 
simply the HJPDE. 



4 Formulation for singular second order Lagrangians 



We consider now the application of the formalism developed in the previous sec- 
tion to a system with a singular second order Lagrangian. When the Hessian 
matrix has a rank n — R, R < n, the momenta variables will not be independent 
variables among themselves. In this case we can choose the order of accelerations 
Qi=(ii in such a way that the minor of rank n — R in the bottom right corner has 
nonvanishing determinant: 



(30) 



det 



d 2 L 
dq a dq b 



det 



d-K b 

dq a 



^ 0; a, b = R + 1, n 
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So we can solve the n — R accelerations q a in terms of the coordinates (q, Q 
the momenta ir a and the unsolved accelerations Q a (a— 1, R) as follow: 

(31) Q a = f a (qi,Qi,n,i 

Since the momenta 7r are functions of the accelerations Qi we can substitute 
the expressions ( |3"1~D and obtain: 

(32) Hi = gi (q t , q u q^j = g { (q u q u f a , q a ^j 



(33) 7Tj — (ji yqii qi, 7r a , q c 

Since we have n a = g a the other n — R functions g a can't contain the unsolved 
accelerations q a or we would be able to solve more of the accelerations as 
functions of the other variables, which contradicts the fact that the rank of the 
Hessian matrix is n — R. So, we can write for the momenta ir a : 

(34) tt q = -HI [q u q { ,p a , 7r a ) 

We can obtain a similar expression for the momenta p a : 

(35) p a = -HP (q h q t ,p a , 7r a ) 

Anyway, from a general constraint given by any expression like equation (§) 
we can always obtain expressions like equations ( |3"1| ) and ( |5oD (see ref.|TT|). 



The Hamiltonian Hq, given by equation fl28[ ), becomes 



(36) H = p a q a + 1 a P a \ pp= _ H P + T^afa 

+ n a \ ni3= _ H K - L \ q h q a , q a = f a 



where a, (3 = 1, R; a — R + 1, n. On the other hand we have 

^ dH df a dL dL df a 
(37) — — = ir a ^— + tt q 7 7 — = 

dq a dq a dq a dq a dq a 
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so the Hamiltonian H does not depend explicitly upon the accelerations Q a . 

From this point we will adopt the following notation: the coordinates t and 
q a will be called t and t a , respectively, and q a will be called i a . The momenta 
p a and 7i a will be called P? and P£, respectively, and the momentum P Q will be 
defined as: 

(38) P = f 

Then, to obtain an extreme of the action integral, we must find a function 
S (qi, <ji, t) that satisfies the following set of HJPDE: 

(39) H' = P + H ( t ,t a ,i a ;q a ,q a ;p a = — ; n a = ) = 

(40) = F» + , 

(41) H': = PI + HI \to, t a , t a ] q a , q a ; Pa 
From the definition above and equation ( |3"6D we have 

07l b Q q a OK b 01T b 

dm _ dm - am 

43 — ^ =q b - q a — q a 

01X b OTT b OTT b 

and 

dH' Q - dm a _ dHg - - - dm a _ dm a 

where a = 1, R; a, b — R+ 1, n. 

Remembering that <Zi=?i and multiplying by dt = dt we have from equations 
flU and (H) 





dS 


dS \ 


q a ;p a -- 


dq a 


~dq a ) 




dS 


dS N 


q a ;p a 


dq a 


dq aJ 




dS 


dS N 


,q a ;p a 


dq a 





, _ dH' dm , , _ 

45 d q b = —^dt + -rr^dt a + —^d t a 

dir b dir b drc b 

dHL , dm , dm 

(46) dq b = —^dt + —^dt a + -^d t a 

dp b dp b dp b 
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and also 



(47) 

(48) 

(49) 
for (3 

(50) 
(51) 



dq a 



dm 
dPl 



\dt Q + 



dm 



dP 



p dtp 



~dPJ 



■d tp~- 



d_Hl 
dPg 



dtr- 



5 a pdtp 



dt r 



dq c 



dm oh' p dm; _ oh'p 

:dt + ^r&dt 3 + ^rd tg= ^r-dt 



dPJl 



BP* 



dPl 



$apd t/3— d t c 



dH' 

dqo = dt = -zr^-dt 



dHf 



-dtg + 



dm 



-d tp~- 



dm 



°dt 



dp dp dp p dp 

.,R; so we can write the equations (EH) and (|46"D as 



, - dH' dH'P , dH'; , _ 
d Qi= ^r^dto + ^^dt a + ^^d t c 



dtr 



dqi 



d-Ki 



dH' , 
°-dt 



d-Ki diii 
dH'P j dH'^ _ . 



1, ...,n 



n 



dp, dp* dpi 

If we consider that we have a solution 5 (<&, t) of the set of HJPDE given by 
equations (]39f) , (ffQ|) and (EJ) then, differentiating those equations with respect 
to 9j, we obtain 



(52) 

(53) 
(54) 



,„ dH' Q d 2 S dH'y d 2 S dH' d 2 S 

Wi + dP dtd + d Pa dq a d % + d7C a dQadQi 

dH'P dH'P d 2 S dH'P d 2 S dH'P d 2 S 







dq q 

dm; 



+ 



dPg dtpd q t 
d 2 S 



dm; 



dp a dq a dq { 

32 i 



dm; d 2 s dm; 



dn a dq a dq~i 
d 2 S 



d Qi ' dP$ dipdQi dp a dq a d q % dn a d % d q % 
whereas for qi we get 



(55) 

(56) 
(57) 



dH' t dH' Q d 2 S dH' Q d 2 S ^ dH' d 2 S 



dqi dP dtdqi 



dH'; dH'P d 2 S dH'; d 2 S 



dp a dq a dqi dr[ a d q a dqi 
dH'P d 2 S 







dqi 

dm; 



+ 



dP| dtgdqi dp a dq a dqi dn a d q a dqi 

dm; d 2 s ^ dm; d 2 s ^ dm; d 2 s 



dqi dPfi d tg dq { dp a dq a dqi dir a d Q a dq { 
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and for t we have: 

dH' Q d% d 2 S dH' Q D 2 S dH' d 2 S 



(59) ggg | dH 2 d " S | ^ p g2g | dH « d " S Q 

<9t dP$ dtpdto dp a dq a dt dn a d q a dt 

(60 ) dH ? | dH ? d " S | ^ d " S | ^ ^ -Q 

dt dPfi d tp dt dp a dq a dt dn a d q a dt 

Making Z = S <ji, t)and using the momenta definitions together with equa- 
tions (^) and fl5T|) we have 

ir7 dS , dS 1 dS T _ 95 , , _ 

(61) dZ = ^7^0 + «T"«*a + o"=^" *a «9a + TrE~ d 9 a , 

ot ot a v ta oq a o q a 



(62) 



dZ 



'dm , OH'p , dH£ 



9p a 



9p a 



5p a 



dH% ^ 

— flto "I - — 7. dt a -\- —- d ta 



(63) 



dZ 



-H + Pa —^ + 7T a 



dp a 

+ {-HI + Pa 



dir a 



dtr 



8H* . 



+ [-H a -+ Pa - 



dp a 

oh': 

dp a 



+ 7Ta + 



d t c 



and, from momenta definitions: 

, x ,„ d 2 S , d 2 S , d 2 S , <9 2 S <9 2 S , _ 

(64) dP = —dt + —dt a + m d qa + ^J-d ta +^-d % 
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, , , d 2 S , d 2 S , d 2 S , a 2 5 ,_ 
(65) apj = — — ^ato + — — —dt a + — — - — dq a + — — 7^-d t a 



d 2 S 



dqidt 



(66) 



d-TTi 



dqidt a 
d 2 S 



dqidq a 



dqid t a 



dqid q, 



-dq a 



d 2 S , d 2 S , 

dto + - , ata + - _ - — dg a 



dq~idt dQi dt 
d 2 S 



+ 



"rf + 



, d Qi dq a 

d 2 S 



-dq a 



d^dta dq t dq a 
Now, multiplying equations (|58|) by (i£ , contracting equations fl59|) and (|60| 
with (it a and d t a (respectively) and adding them all to equation ( pi] ) we get: 



(67) dP Q + ^dt + ^^dt a + ^^di 



dtn 



d 2 S 
dt dq a 

d 2 S 
dt d q a 



dt dt 

am , dH'i 

dq a — dt 



dp a 



. . dt a 

OPa 



d_K 

dp a 



' dH' dH'i , 
d q a -^-dt - ^r^dt r 



dir a 



^di 

07T a 



In a similar way, using the same steps with equations (|55| - |57|) and (|65D we 
obtain: 



(68) 



dft + -TT A dt + —^dt a + —^dt 

oqi dqt dqi 

d 2 S 



dqidq a 
d 2 S 
dqid q a 



dq a - 



am 



dp t 



'-dtr 



dH'P 

dp a 



' , - dH' 
d q a —z-^dta 

OTX a 



dH' p 
a dt 

dHa 



d_K 

dp a 

dH'i 
dn a 



d t c 



d 



And, finally, using equations (|52| - ^4[) and (|66|) we have 



dH' dH'P , dH'I , - 
69 diii + — ^dt + — -?-dt a + — S-d t a ~- 

dq { dq { dq { 
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d 2 S 
d Qi dq a 

d 2 S 
dQidq* 



dq a — - — aU 



dq a 



dp a 

dH^ 

dlXa 



dto - 



dp a 



dp a 



°Ldt - d -^-dt 
OTl n 07T ri 



If the total differential equations given by (f45j ) and (fig) are valid the equations 



fl68|) and (|69|) become: 



(70) 

(71) 
(72) 



dP n 



dpi 



dm 



°dtn 



dH ~ P -dt„- d Sdi c 



8tr 



dt 



dH' , OH'p , dH'* , _ 



dqi 



dlXi TT^r-dtv, 



dqi 

dm 



dt r 



dqi 
dH'^ 
d^ 



d t c 



dQi d^ 

These equations together with equations (fE?D, Q5DD, fl5"T| ) and ( |B3| ) are the total 
differential equations for the characteristics curves and, if they form a completely 
integrable set, their simultaneous solutions determine S ((ft, to) uniquely by 
the initial conditions. 



5 Example: Podolsky generalized electrodyna- 
mics 

In this section we will consider a continuous system with Lagrangian density 
dependent on the dynamical field variables and its derivatives upon second or- 
der: C = C (■0, dip, d 2 %fj). We adopt the metric rj^ = diag(+l, —1, — 1, —1) with 
Greek indices running from to 3 while Latin indices run from 1 to 3. As stated 
previously, the generalization of the formalism presented in Sections 4 and 5 is 
straightforward, being necessary only to consider that the Euler-Lagrange equa- 
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tions of motion are now given by 

(73) w~ d * {w^j) + d " du = 

■a ■■ a 

and that the momenta, conjugated respectively to ip and ip , are: 

dC n ( dC \ ( dC \ 

(74) V a = —a ~ 2d k \ —. raT ~ <% 



<9 *0 "0 )/ \8 if) / 

(75) vr a = — — 

d if) 

The Hessian matrix is now: 

d 2 C 

(76) H ab = —^— b 

d if) d if) 

With these modifications we consider the case of Podolsky electrodynamics 
which is based on the following Lagrangian 

(77) C = -\f^ + a 2 d x F aX d p F a ? 

were F^ = d ll A v - d u A^. 

An analysis of the Hamiltonian formalism for this theory was carried out in 



ref. ||12|| and we compare some of the results presented there with the formalism 



developed here. The Euler-Lagrange equations are 

(78) (l + 2a 2 n) d p FJ = 

with our dynamical variables chosen as A^ and ~fit=)t ■ The conjugated momenta 
given by definitions ([74]) and ([75]) are: 

(79) Vii = -F 0ll - 2a 2 (d k d x F 0X 5 k „ - d d x F„ x 

(80) Ti ll = 2a 2 (d x F 0X 5\-d x F ll 
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The primary constraints are: 



$ x = 7T « 



32) $ 2 = po - d fc 7r fc « 

_ t 

Using the definition of 7r we can write the accelerations A as: 
33) 

The canonical Hamiltonian is given by: 



.A = ^ + d k F tk + d l Ao 



34) 



H c = / d 3 x 



Vp, A 1 * +7r M A -C 



Using equation fl83|) we get 



35) # c 



ffx 



A° d l 7r { + Pl A 1 + -^W + TCidkF ik + Ti.d 1 Ao +]f imu F^ 
4or 4 



+- (A -Mo) (I' -dM ) - a 2 (d k A k -d k d k A ) (d { t -d i d i A i 



According to Dirac's formalism the total Hamiltonian is: 
36) H T = H C + J d 3 x {C x (x) $! + C 2 (x) $ 2 ) 

The consistency conditions result in: 



37) 



$1 = {$ 1; if T }«0 



$ 2 ={$ 2 ,# T } = aV/c~o 

So we have a secondary constraint given by 



39) 



$3 = dVfc ~ 
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and the consistency condition results in $ 3 = {$ 3 ,i^} w 0. All constraints are 
first class so the extended Hamiltonian is: 

(90) H E = H C + J d 3 x {C x (x) $i + C 2 (x) $ 2 + C 3 (x) $ 3 ) 

The equations of motion for the dynamical variables, given by A = {A a , He}, 

are: 

(91) A°=A° +C 2 - A=A l -d l C 3 

This simply means that A is defined as A plus additive arbitrary functions. 
Besides, A = [A ,H E j gives 

(92) A°= C i; A = ^ + d k F' lk + &> Aa, 

which mean that both A° and A are arbitrary while we obtained again equation 



31). 

For the momenta variables 7Ti— {ni, He} and P a = {p a , He} give: 

(93) 7r i= -F 0i - 2a%d k F k - Pi 

(94) P =-d i F 0i -2a 2 d i d i d k F k 

(95) V l = -did k 7r k + d k d k 7Ti - d k F k t 

Equation ( p3|) is the definition of Pi given by equation (|79|) and together with 
fl9~j| ) it gives constraint $ 3 . 

Now, using Hamilton- Jacobi formalism we have: 

dS 

(96) H' = H C + P ;P = — 

(97) ifj = vr ; H' 2 = p - d k ir k 
The total differential equation for A* is 

(98) dA* = 9 -^dt + 9 -^d Ao +^dA 

dpi dpi dpi 
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(99) 



dA = -z-±dt = -z-^dt => dA =A l dt 



dpi dpi 

which is completely equivalent to equation fl9~T|) since C3 is arbitrary. For A* we 
have: 
(100) 

un i un j 

(± 



, -i dm , dm - am 1A am , oh c , 

d A = -rr^dt + —±d A +^dA = —^dt = —±dt 

UHj OTXj OTtj OTTi OTXj 



(101) 



<l A-- ^rr l + d k F lk + d l A ) dt 



Again we have a result in agreement with Dirac's method result given in fl92| 
For the momenta pi and po we have 

, inn , , , am , aH[ , - am 1A am , aH c , 

^ dp = ~id dt - ai d A ° ~id dA « = -ai dt = ~id dt 



(103) 

(104) 

and 
(105) 



dp % = - / d 3 x 



' dP k \ _ 1 in dF jn 
dAi i 2 dAi 



dt 



dp 1 = -a i d k n k + d fc av - a k F 



dt 



dp° 



am , 9m , _ 

-dt - — — ~d Ao 



dA 



dA d A, 



dH *dA = A 



dA 



(106) dp 



d x 



A -a' a, 



a (a -diA 



-2a 2 (di A* -di&Ao) 



dA 

d (d k A k -d k d k A 
dA 



dt 



(107) dp = [-di (X -d'Ao) - 2a 2 d k d k (di A 1 -did'Ao 

Finally for n % we have: 



dt 



dHL , dH[ , - dH!) , , 
108) aV = -—^dt - —^-d A -—^dA c 



d Ai d Ai 



dAi 



dH '»-dt = - d -^dt 



dA 



dA 
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(109) 



dn l 



- / (Fx 



P 



idAi 



d Ai 

2a 2 (d,j A j -dj&A 



+ (X j -&Ai 



d (Aj -djA 
dAo 



d (d k A -d k d k A 
dAo 



dt 



(110) 



-p l - F 0i - 2a z d l d k F 



2»a J?0k 



dt 



Equations ( p^D , ( |107| ) and (|TT0|) are completely equivalent to (|9|), (||) and 
(B5J) ; consequently equations ( |107| ) and ( |110| ) give us the secondary constraint 
that isn't present in the total differential equations. 



6 Conclusions 

We obtained a generalization of Hamilton- Jacobi formalism whose results agree 
with those obtained using Dirac's formalism. In this formalism those coordinates 
whose correspondent accelerations can't be solved in function of the momenta 
are arbitrary variables of the theory. We obtained a set of Hamilton- Jacobi 
partial differential equations in terms of these variables and from this set we 
obtained the equations of motion of the system as the total differential equations 
for the characteristics. These total differential equations so obtained must satisfy 
integrability conditions and for these conditions to be satisfied the nature of the 
constraints (first class or second class) will play an essential role. The study of 
these integrability conditions is in progress as well as the generalization of the 
present formalism for Lagrangians of order higher than two. 
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